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MAGNETOTHERMAL OSCILLATIONS IN PRESSURE- 
ANNEALED, PYROLYTIC GRAPHITE 
by Dennis J. Flood 
Lewis Research Center 

SUMMARY 

Magnetothermal oscillations have been observed for the first time in a sample of 
pressure-annealed pyrolytic graphite. Changes in sample temperature of the order of 
5x10 K were measured using a field modulation technique. The period was studied as 
a function of the angle between the c-axis of the crystal and the magnetic field direction. 
Two periods were observed whose values for magnetic-field strength H parallel to the 

_ i 

c-axis are 0. 151±0. 005 and 0. 85±0. 20 tesla . The former agrees well with the period 
of the majority hole carriers obtained by other workers using different techniques in both 
single crystal and pyrolytic graphite. The second period is tentatively ascribed to the 
minority electrons and is lower than the values previously reported 

The measurements reported here have extended the range of magnetic field over 
which quantum oscillations in any equilibrium phenomenon have been observed in pyro- 
lytic graphite. The phase factor y has been determined, with the use of the high mag- 
netic field data, to be 0. 50±0. 10. 


INTRODUCTION 

The term " magneto thermal oscillation" refers to the oscillatory dependence of the 
temperature of a thermally isolated sample on inverse magnetic field. Kunzler, Hsu, 
and Boyle (ref. 1) originally observed the phenomenon in a single crystal of bismuth by 
a simple dc technique. Following that, LePage, Garber, and Blatt (ref. 2), using a vari- 
ation of the dc technique of Kunzler et al. , reported the effect in beryllium. McComb and 
Seidel (refs. 3 to 5), using a field modulation technique, investigated the effect in bis- 
muth, antimony, and beryllium. 

This report discusses the observation of magnetothermal oscillations in pressure- 
annealed, pyrolytic graphite by means of both a field modulation technique and the modi- 


fied dc technique of LePage et al. The oscillations were studied as a function of orienta- 
tion in fields ranging from 0. 35 tesla to 10 tesla at temperatures between 1. 3 and 1. 5 K. 

Pyrolytic graphite is essentially a polycrystalline material which exhibits a high 
degree of order in one direction, and nearly complete disorder in the plane perpendicular 
to that direction. The behavior of the electronic properties of pyrolytic graphite can 
qualitatively approach that of the electronic properties of single crystal graphite as the 
order is improved through the high-pressure annealing process (refs. 6 and 7). 
Williamson, Foner, and Dresselhaus (ref. 8) have observed deHaas-van Alphen oscilla- 
tions in pyrolytic graphite, and compared the results with those obtained from deHaas- 
van Alphen measurements in single crystal graphite by Soule (ref. 9) and with magneto- 
resistance measurements (Shubnikov - de Haas effect) in single crystal graphite by Soule, 
McClure, and Smith (ref. 10). They found that the general features of the majority 
carrier pockets are identical for single crystal and pyrolytic graphite, but that details 
for the minority carrier pockets near the zone corners differ somewhat. Their conclu- 
sion was that the Slonczewski- Weiss band model for graphite (ref. 11) is applicable to 
pyrolytic graphite, and that the electron and hole Fermi surfaces are closed. 


THEORETICAL 

Magnetothermal Oscillations 

According to thermodynamics, the change in internal energy U of a homogeneous 
system in which temperature T, magnetization M, and magnetic field strength H are 
the only state variables of interest is given by 


dU = d'Q - d'W = C H dT + T dH 


3T / 


(1) 


H 


where d'Q is positive for heat flow into the system, d'W is positive for work done by 
the system, and Cjj is the total heat capacity at constant magnetic field. The rate at 
which heat flows from the system into a reservoir is 


d’Q 

dt 


-/c(T - T ) = -k AT 
v s r 7 


( 2 ) 


where T is the system temperature, T is the temperature of the reservoir, and k 
s r 

is simply Ka/L, where K is the thermal conductivity of the path of length L and cross- 
sectional area a which connects the system to the reservoir. From equations (1) 
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and (2), the temperature difference between the sample and the reservoir AT, is given 

by 


- K AT 


dT 


'H 


dt 



(3) 


The magnetization of a system of conduction electrons contained within a solid has 
been derived by Lifshitz and Kosevich (ref. 12). Their original calculation has since 
been refined by several authors (refs. 13 to 15). The oscillatory part of M, which is 
all that is of interest here, is given by 



where 


= " 2Vk B 



9 2 A(E f ) 


3k„ 


- 1/2 
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ni gm* 
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27T m*ck 
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Dingle temperature: 


de Haas - van Alphen frequency: 


x = 


K 

ffk B r D 


f =£ A o<V 


(5c) 


(5c) 


(5d) 


k Q defines the location in k-space of the extremal area denoted by A Q (E F ), E F is the 
Fermi energy, m* is the effective mass given by (h 2 /27r) 9 A(E f )/3E f , V is the volume 
of the sample, g is the spin- splitting factor (ref. 13), r D (ref. 14) is the average time 
an electron travels between collisions, and kg is the Boltzmann constant. The other 
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symbols are defined in the appendix. The quantity y which appears in the argument of 
the sine in equation (4) is the phase factor which appears in the expression developed by 
Onsager (ref. 16) for the quantization of the area of the orbits traveled by the electrons: 

A =( n+ y)®M (6) 

he 

where n is an integer quantum number which labels the energy levels (Landau levels) 
of the electrons in the magnetic field. For free electrons, or for electrons which obey 
a general quadratic dispersion law, y is equal to 1/2. (This is true basically because 
the Schroedinger equation in such instances reduces to that for the harmonic oscillator, 
with eigenvalues E n = (n + l/2)Ho>*). Chambers (ref. 17) has recently pointed out that 
when interband transitions can be ignored, that is, only weak fields and weak electron- 
electron interactions allowed, y is expected to be 1/2 even for arbitrary Fermi surfaces. 

The expression in equation (4) is subject to some limitations. The magnetic field 
and temperature must be such that Ep »kgT and (heH/m*c) »kgT. It is also 
necessary that the argument of the sine term in equation (4) be large, that is, that 
(2rf/H) » 1. The last condition means simply that the expression for M is no longer 
valid in the so-called "quantum limit regime, ” where the Landau-level spacing becomes 
comparable to the Fermi energy Ep. 

Direct substitution of equation (4) into equation (3) results in 
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The expected amplitude of the oscillations is at most of the order of 10 K, and for 
many materials the amplitude will be less than 10 -4 K. For near-adiabatic conditions, 
the term on the left in equation (7) can be neglected, and the change in sample tempera- 
ture is given by 
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The last term on the right in equation (8) includes all the dissipative effects acting to 
change the temperature of the system such as joule heating in the thermometers, eddy 
currents, and so on. 


Band Theory of Graphite 

The crystal structure of graphite is well known and is illustrated in figure 1. The 
lattice is hexagonal close-packed, and the hexagonal Brillouin zone corresponding to it is 
shown in figure 2. Only a small fraction of the zone is occupied by the current carriers, 
so that only the energy bands near the vertical edges of the zone are important for the de- 
termination of transport effects. Accordingly, Slonczewski and Weiss (ref. 11) have de- 
rived an effective-mass Hamiltonian for the energy bands near the vertical edges. Their 
original calculations were developed further by McClure (ref. 18), who included the effects 



Figure 1. - Crystal structure of graphite. Centers of hexagons in one layer lie over corners 
of hexagons in layer beneath it. Unit cell parameters are a Q = 2. 456 A and c 0 = 6. 712 A. 
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Figure 2. - Brillouin zone of graphite (from Williamson et al. (ref. 15)). 

of a magnetic field. The energy bands derived by the former are shown in figure 3, which 

indicates a hole pocket centered at k = 0, and two electron pockets with extremal areas 

near c k /2 tt = ±0. 35, where c is twice the spacing between layers of the crystal, 
o z o 

These electron and hole majority pockets arise from the energy band labeled Eg. In ad- 



Figure 3. - Energy bands of graphite near vertical zone edge HKH at zero field as 
given by Slonczewski-Weiss model. Small pocket of electrons shown at right 
side of figure arises from Ei band and joins smoothly with pocket of electrons 
formed by Ej band (from Williamson et al. (ref. 15)). (Small pocket forms caps 
of ellipsoids at points H in fig. 2.) 
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dition there are small minority pockets of electrons near the points c 0 k z /27r = ±0. 5 which 
arise from the Ej band. The minority pockets join smoothly with the electron pockets 
in the Eg band in the absence of a magnetic field, or when spin-orbit interactions are 
neglected. Figure 4(a) shows the Fermi surface in a reduced zone scheme for the case 




(a) Minority electron pocket at 
zero field, neglecting spin- 
orbit interaction. The orbit 
at k z « rrc 0 is not stationary. 


(b) Effects of spin-orbit inter- 
action. Two stationary orbits 
are produced as a result of 
the separation of the Fermi 
surfaces of the two pockets. 


Figure 4. - Details of Fermi surface near the Brillouin zone corners in the extended zone scheme 
(from Williamson etal. ref. 15). 


where the two pockets are joined smoothly onto one another, and figure 4(b) shows the 
effects of including a spin-orbit interaction in the Hamiltonian. Only for the second case 
does a true extremal deHaas - van Alphen orbit exist. (McClure (ref. 18), however, has 
shown that the magnetic energy levels, even in the absence of a spin-orbit interaction, 
exhibit stationary points at the same values of k z as in the preceeding discussion.) This 
picture may not be complete. Soule (ref. 9) has suggested, on the basis of acceptor- 
impurity doping studies in single crystal graphite that there may be some minority hole 
pockets as well. These may be located around the majority hole pocket as shown in fig- 
ure 5. 
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Figure 5. - Brillouin zone with possible location of minority hole 
"outrigger" Fermi surfaces as reduced about zone edge (from 
Soule, ref. 9). 


EXPERIMENTAL 

Experimental Techniques 

There are basically two ways to detect magnetothermal oscillations: resistance 
bridge thermometry (refs. 1 and 2), and afield modulation technique (refs. 3 to 5). Fig- 
ure 6 contains a schematic diagram of the experimental setup for the former. A low- 
frequency ac resistance bridge is used to detect the change in resistance of a carbon 
thermometer attached to the sample. The bridge input voltage is supplied by the refer- 


Resistance 

bridge 


Lock-in 

amplifier 


X-Y recorder 










Figure 6. - Schematic diagram of bridge technique. 
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ence oscillator of a lock-in amplifier. The output of the bridge is then synchronously de- 
tected by the lock-in amplifier and converted to a dc voltage which is fed to an X-Y re- 
corder. This technique requires a very careful adjustment of heat leaks from the sample 
to the bath so that eddy current heating in the sample will not raise its average tempera- 
ture appreciably. It is also necessary, however, to provide a long enough thermal time 
constant so that the sample can vary its temperature at least enough to be observed. 
Thermometer magnetoresistance can also affect the measurements in the form of a large 
monotonic background. Nonetheless, in particular situations the technique is quite useful 
and can produce calibrated amplitudes for the oscillatory temperature changes. 

A schematic diagram for the field modulation technique is shown in figure 7. A 
small- amplitude, low-frequency, modulating magnetic field is superimposed on the mono- 
tonically swept dc field. The resulting ac component of the temperature change is mea- 
sured by a lock-in amplifier. In this case the time constant from sample to bath must 
only be long compared with the period of the modulating field, or whatever harmonic of it 
which is being detected, and can be made short enough so that the sample remains at an 
average temperature very close to that of the bath. It is necessary, however, for the 
thermal relaxation time between thermometer and sample to be short compared with the 
period of the modulating signal so that the thermometer will follow the sample tempera- 
ture at the detection frequency. This is usually not a difficult requirement to satisfy. 
Changes in average sample temperature which arise from eddy currents induced by the 
slowly swept dc field will therefore be greatly reduced under these conditions. Moreover, 
eddy currents induced by the modulating field will cause heating effects primarily at 
twice the frequency of the ac field and are eliminated by detection at the fundamental. 
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Figure 7. - Schematic diagram of field modulation technique. 
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Sample and Apparatus 


The sample was a rectangular rod 3. 18 by 3. 18 by 12 millimeters (1/8 by 1/8 by 1/2 
in. ) of pyrolytic graphite and was cut with its long dimension perpendicular to the c-axis. 

It had a resistance ratio R^qq/R^ g °f 1/2 for current flow in the basal plane. 

The sample holder is shown in figure 8. The spiral gear was machined such that a 
complete rotation by it of 360° caused the gear on the sample bed to advance by one tooth, 
which in this case corresponded to 10°. The drive shaft on the spiral gear was a length 
of stainless -steel tubing which was sealed at the top end, fed through the vacuum line, 
and bolted to a brass rod at the bottom. The brass rod passed through the top of the sam- 
ple chamber, which acted as a heat sink. The spiral gear was fastened in place at the 
bottom of the brass rod with a setscrew, and its position could be adjusted to obtain a min- 
imum amount of backlash. The apparatus provides one degree of freedom in an external 
magnetic field. 

The open table on which the sample was mounted was machined from oxygen-free 
high- conductivity copper and was supported by bushings of the same material which were 
soldered into the brass wall of the chamber. Wires to the sample could be clamped to 
the edges of the table so that rotating the holder would not cause the leads to pull directly 
on the sample. The clamp also provided a second heat sink for the leads in addition to one 
located on the wall of the chamber. The sample was secured to the table with magnet- 
wire varnish and a second clamp. The thermometer was glued on the end of the sample 
away from the clamp with varnish and was electrically insulated from the sample by a 
piece of cigarette paper. 

The dc magnetic field was provided by a 6. 35-centimeter (2-^-in. ) bore, 10. 5-tesla 
superconducting solenoid. The field was measured by a copper magnetoresistor which 
had been calibrated by nuclear magnetic resonance at several values of field strength. 

The estimated error in the calibration is ±0.2 percent. The modulating field was provided 
by a coil wrapped around the outside of the vacuum chamber containing the sample. Mod- 
ulating fields up to about 0. 01 tesla were obtainable. 

An insert Dewar with a 5. 33-centimeter (2. 1-in. ) inner-diameter tail was used to 
reach temperatures below 4. 2 K. Temperatures near 1 K could be attained when no power 
was supplied to the modulation coils. Measurements were typically made in the vicinity 
of 1. 3 K. 


Experimental Procedure 

Although both methods were tried, the field modulation technique was used almost 
exclusively. The primary difficulty with the bridge technique is the large monotonic back- 
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Figure 8. - Sample holder detail. 
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ground arising from the thermometer magnetoresistance. This can be reduced by com- 
paring two resistors with each other, but the necessary compensation is difficult to 
achieve because of differing field and temperature characteristics . This problem is not 
encountered with field modulation, since to a first approximation 

AR /aR\ /3T\ JdR\ (9 

AH \3T/ h \aH/ g W T 

The second term is very nearly constant over the entire field range and provides only an 
offset from zero, not an increasing background. Figure 9 is a plot of resistance as a 



function of magnetic field strength under isothermal conditions and shows that 
(SR/SH)^ ^ is relatively small and does not change much over the field range used. 

The variation of R with temperature for the same 100-ohm (nominal room temperature 
value), 1/8-watt resistor is shown in figure 10. In the region around 1. 3 K, (3 R/ST)jj 
is very large, and for small changes in temperature, nearly constant. Hence, changes in 
AR/AH are the result almost entirely of changes in (9T/3H) . 

O 

The thermal response of the system was such that temperature changes could be de- 
tected at frequencies up to 33 hertz. Amplitudes at the higher frequencies were reduced 
considerably. Typically, measurements were made at 4 hertz. The integrator time con- 


12 




stant at the output of the lock-in amplifier was usually 3 seconds. (As mentioned earlier, 
the frequency response is limited by the thermometer-sample contact. After the present 
data were collected, attempts to improve the response were partly successful, enabling 
data to be taken at 33 Hz with almost the same sensitivity as before. ) 

The oscillations were studied as a function of angle between H and the c-axis of the 
crystal by manually rotating the sample holder to the desired position and then sweeping 
the magnetic field. The angular position of the sample was recorded by measuring the 
voltage across a 10-turn potentiometer. The initial error in position was estimated to be 
±1°, and the relative error of any setting was ±0. 5 percent. 

A dc current of the order of 1 microampere was passed through the thermometer, 
and the voltage changes measured were in the range from 1 to 20 microvolts. Thermo- 
meter sensitivity was typically 0. 2 volt per degree kelvin for the operating points chosen 
(i.e. , T “ 1. 3 - 1. 5 K). 


Results and Discussion 

The angular variation of the period of the observed oscillations is shown in figure 11. 
The value for H parallel to the c-axis is 0. 151±0.005 tesla - '*'. This is in good agree- 
ment with the hole period obtained by Soule et al. (ref. 10) from Shubnikov - de Haas 
measurements in single crystal graphite, and with that obtained by Williamson et al. 
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■ -20 0 20 40 60 80 

Angle between H and c-axis, deg 

Figure 11. - Variation of period with angle. 

CS -49783 

(ref. 8) from de Haas - van Alphen measurements in pyrolytic graphite. The uncertainty 
in the period reported here is a direct result of the difficulty in determining the exact 
location of the small, lower field peaks. The field at which any peak occurred could be 
determined from the X-Y recorder to within ±0. 015 tesla. The period as a function of 
angle was determined by rotating the sample to the new position and sweeping the field. 
Consequently, the same uncertainty in period is associated with each point in the plot of 
figure 11. 

The hole oscillations have superimposed on them a slowly varying envelope curve 
(fig. 12). The 1/H plot in figure 13 gives a period of 0. 85±0. 20 tesla” * . This os cilia - 



Magnetic field strength, H, T cs-49784 


Figure 12. - Recorder trace of magnetothermal oscillations showing 
superimposed envelope curve and estimated locations of maximums and 
minimums. 0= 0° from c-axis; temperature, 1.32K. 


14 



Figure 13. - Peak number as function of lfH for the long-period en- 
velope oscillation. Period, 0.85±0. 2 tesla' 1 . 


tion is tentatively ascribed to the small pocket of minority electrons located at the points 
H, H' in the Brillouin zone (see fig. 2). The assignment is only a suggestion. William- 
son et al. (ref. 8) found a minority period of 1. 35±0.01 tesla -1 in single crystal graphite, 
and 2. 24±0. 03 tesla - 1 in pyrolytic graphite, for the same orientation. On the other hand, 
Soule (ref. 9) has reported the possible existance of three minority hole pockets in single 
crystal graphite which may be located around the majority hole pocket in "outrigger" 
fashion (see fig. 5). The period of this minority carrier oscillation varied from 0. 58±0. 1 
to 1. 35±0. 3 tesla -1 , depending on the amount of acceptor impurities introduced into the 
single crystal samples. The period of the superimposed oscillation given here could not 
be determined with very high precision. The presence of a remnant field in the super- 
conducting solenoid prevented acquisition of data below about 0. 35 tesla. Hence, only a 
few maximums and minimums could be observed, and the superposition of the shorter 
period hole oscillations made their location difficult to determine. The envelope oscilla- 
tion was observable only at the lowest temperatures, and an attempt to investigate the low 
field region in another magnet failed because the same temperatures could not be reached. 

The question of whether the observed long-period oscillation comes from a hole sec- 
tion or an electron section cannot be answered here. The sample has not yet been anal- 
yzed for impurities, since it was felt advisable to complete other measurements of its 
transport properties before cutting it. The presence of acceptor impurities would lower 
the period for a minority hole pocket, while the presence of donor impurities would lower 
the period for a minority electron pocket. A systematic investigation of the effects of 
such impurities would be necessary to establish the type of carrier. A spectrographic 
analysis of a similar piece of pyrolytic graphite has indicated the presence of 30 ppm of 
aluminum, and only trace amounts (<1 ppm) of boron. Aluminum would be a donor im- 
purity and boron an acceptor. The presence of the aluminum favors the assignment of the 
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oscillation to the minority electron pocket, since the samples were presumably prepared 
under the same conditions. The majority volumes are relatively unaffected by small im- 
purity concentrations. 

Figure 14 shows a recorder trace of dT/dH as a function of H up to 8 tesla. The 


dV/dH 


l l i i i [ i i i 

0 1 2 3 4 5 6 7 8 

Magnetic field strength, H, T 

Figure 14. - Recorder trace of magnetothermal oscillations. 0 = 17° from c-axis ; tempera- 
ture, 1.4 K. 

magnetic field is tilted 17° from the c-axis. The figure clearly indicates how quickly the 
low quantum number regime is reached for the hole oscillations. The waveform is begin- 
ning to show some distortion by about 1. 0 tesla, and is definitely nonsinusoidal for fields 
above 2 tesla. The frequency of the oscillations in figure 14 is 6. 95±0. 20 tesla. Since 
the peaks in the figure actually correspond to maximums in dV/dH, where V is the 
voltage across the thermometer, they are shifted from the maximums in dT/dH by 180°. 
(Resistance as a function of temperature has a negative slope.) The condition for the 
occurrence of the n u maximum is therefore 



£11 _ 2iry - - = 27m - - (10) 

H 4 2 

Table I lists the peak positions predicted by equation (10) and the fields at which they were 
actually observed. Equation (10) is not expected to be valid in the so-called quantum 
limit. The agreement between observed and predicted maximums is within experimental 
error for fields up to 2 tesla. 

A plot of n as a function of 1/H for the field parallel to the c-axis is shown in fig- 
ure 15. The integers were assigned to the peaks on the basis of the preceding discussion. 
The intercept may be used to calculate the quantity y in equation (10). From the figure, 
the intercept of the line is n = -0. 37±0. 10. Inserting this in equation (10) gives 
y = 0. 5±0. 1. 
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TABLE I. - COMPARISON OF PREDICTED FIELD 


POSITIONS OF MAXIMUMS IN dV/dH 


WITH OBSERVED FIELD POSITIONS 


Peak number, 

Field position, T 


Predicted 

Observed 

1 

5.06 

7 . 8 ± 0 . 02 

2 

2.93 

3 . 2 ± 0. 02 

3 

2.06 

2 . 05 ± 0. 02 

4 

1 . 59 

1 . 54 ± 0 . 02 

5 

1.29 

1 . 27 ± 0 . 02 

6 

1.09 

1 . 08 ± 0 . 02 

7 

9.43 

9 . 30 ± 0. 02 

8 

8 . 30 

8 . 20 ± 0. 02 

9 

7.42 

7 . 50 i 0. 02 

10 

6.70 

6 . 60 ± 0. 02 

11 

6 . 12 

6 . 10 ± 0. 02 

12 

5.62 

5 . 60 ± 0 . 02 



Figure 15. - Peak number as function of 1/H for hole oscillations 
with H along c-axis. 



CONCLUSIONS 


The periodicity in 1/H of the oscillatory variations in temperature of the sample of 
pressure-annealed, pyrolytic graphite used in these experiments agrees well with that 
obtained for hole carriers by other workers in both single crystal and pyrolytic graphite. 
The measured value of the hole period for H parallel to the c-axis was 0. 151±0. 005 
tesla” From observations into the quantum limit, the phase factor y was found to be 
0. 5±0. 10. 

A long-period oscillation was also observed which is tentatively ascribed to the min- 
ority electron pocket occurring at the corners of the Brillouin zone. The period of this 
oscillation was estimated to be 0. 85±0. 20 tesla” * and is not in good agreement with other 
published values for these carriers. One explanation for the different minority period 
observed in this study is the possible presence of donor impurities in the sample. An 
analysis of a similar piece of pyrolytic graphite has indicated the presence of 30 ppm of 
aluminum and tends to corroborate this suspicion. Small traces of impurities could 
greatly affect the dimensions of the small pockets of minority carriers, while hardly af- 
fecting the majority carrier volumes. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 5, 1969, 

129-02. 
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APPENDIX - SYMBOLS 


A 

A(E f ) 

A 

n 

A o(E f ) 


a 


C 


H 


c 


c 


o 


E 


F 


E 


n 


E 

E 

E 

f 


1 

2 

3 


H 

H, H' 

h,K 

K 

K, K' 

k B 

k o 


L 

M 

M 


osc 


symmetry point in hexagonal face of Brillouin zone 
cross-sectional area of Fermi surface 

area enclosed by orbit of electrons in state labeled by quantum number n 

extremal cross-sectional area of Fermi surface 

cross-sectional area of thermal path from sample to heat reservoir 

specific heat at constant magnetic field 

speed of light in vacuum 

twice the interplanar separation along the c-axis of graphite 
Fermi energy 

energy level of free electrons in a magnetic field 
electronic energy band in graphite 
electronic energy band in graphite 
electronic energy band in graphite 

de Haas - van Alphen frequency of oscillations periodic in 1/H 
effective spin- splitting factor for electrons in a magnetic field 
magnetic field strength 

symmetry points along zone edges of hexagonal Brillouin zone 
Planck’s constant (divided by 27 t) 

thermal conductivity of thermal path from sample to heat reservoir 
symmetry points at center of zone edges of hexagonal Brillouin zone 
Boltzmann constant 

value of k_ at extremal area of Fermi surface 
z 

x-component of wave vector of electron 
y-component of wave vector of electron 
z-component of wave vector of electron 
length of thermal path from sample to heat reservoir 
magnetization of a system of conduction electrons 

oscillatory component of magnetization of system of conduction electrons 



m free electron mass 

o 

m* cyclotron effective mass of electron in a magnetic field 

n quantum number labeling the energy level of an electron in a magnetic field 

Q heat adsorbed or released by a thermodynamic system 

R^qq sample resistance at room temperature 

R 4 g sample resistance in liquid helium 

S entropy of thermodynamic system 

T absolute temperature of thermodynamic system 

T H "magnetic temperature" defined by eq. (5b) 

U internal energy of thermodynamic system 

V volume of system of conduction electrons 

W work done by or on a thermodynamic system 

X Dingle temperature 

y phase factor in energy dispersion relation of electrons in a magnetic field 

r symmetry point at center of hexagonal Brillouin zone 

k heat per second per degree transported along thermal path from sample to 

reservoir 

Tjy average between collisions for an electron in a solid 
w* cyclotron frequency of electrons in a magnetic field 
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